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Low energy heavy-ion fusion reactions are governed by quantum tunneling through the 
Coulomb barrier formed by a strong cancellation of the repulsive Coulomb force with the 
attractive nuclear interaction between the colliding nuclei. Extensive experimental as well as 
■ theoretical studies have revealed that fusion reactions are strongly influenced by couplings 

' of the relative motion of the colliding nuclei to several nuclear intrinsic motions. Heavy-ion 

' subbarrier fusion reactions thus provide a good opportunity to address a general problem 

on quantum tunneling in the presence of couplings, which has been a popular subject in 
the past decades in many branches of physics and chemistry. Here we review theoretical 
aspects of heavy-ion subbarrier fusion reactions from the view point of quantum tunneling in 
systems with many degrees of freedom. Particular emphases are put on the coupled-channels 
approach to fusion reactions, and the barrier distribution representation for multi-channel 
penetrability. We also discuss an application of the barrier distribution method to elucidation 
of the mechanism of dissociative adsorption of H2 molecules in surface science. 
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§1. Introduction 



Quantum mechanics is indispensable in understanding microscopic systems such 
' as atoms, molecules, and atomic nuclei. One of its fundamental aspects is a quantum 

tunneling, where a particle penetrates into a classically forbidden region. This is a 
wave phenomenon and is frequently encountered in diverse processes in physics and 
\Q 1 chemistry. 

. The importance of quantum tunneling has been recognized from the birth of 

quantum mechanics. For instance, it was as early as 1928 when Gamow, and inde- 
pendently Gurney and Condon, applied quantum tunneling to a decays of atomic 
nuclei and successfully explained the systematics of the experimental half-lives of 
radioactive nuclei. 1 )' 2 ) 

^ ■ In many applications of quantum tunneling, one only considers penetration of 

a one-dimensional potential barrier, or a barrier with a single variable. In general, 
however, a particle which penetrates a potential barrier is never isolated but inter- 
acts with its surroundings or environments, resulting in modification in its behavior. 
Moreover, when the particle is a composite particle, the quantum tunneling has to 
be discussed from a many-particle point of view. Quantum tunneling therefore in- 
evitably takes place in reality in a multi-dimensional space. Such problem was first 
addressed by Kapur and Peierls in 1937. ^ Their theory has further been devel- 
oped by e.g., Banks, Bender, and Wu, 4 ) Gervais and Sakha, 5 - 1 Brink, Nemes, and 
Vautherin, 6 ) Schmid, 7 ) and Takada and Nakamura. 8 -* 

When the quantum tunneling occurs in a complex system, such as the trapped 
flux in a superconducting quantum interference devices (SQUID) ring, 9 ) the tunnel- 
ing variable couples to a large number of other degrees of freedom. In such systems, 
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Fig. 1. The internucleus potential between le O and 144 Sm nuclei as a function of the relative dis- 
tance. The dotted and the dashed lines are the Coulomb and the nuclear potentials, respectively, 
while the solid line denotes the total potential. Vb and Rt are the height and the position of 
the Coulomb barrier, respectively. -Rtouch is the touching radius at which the projectile and the 
target nuclei start overlapping significantly with each other. 



the environmental degrees of freedom more or less reveal a dissipative character. 
Quantum tunneling under the influence of dissipative environments plays an impor- 
tant role and is a fundamental problem in many fields of physics and chemistry. 
This problem has been studied in detail by Caldeira and Leggett. 10 ^ This seminal 
work has stimulated lots of experimental and theoretical works, and has made quan- 
tum tunneling in systems with many degrees of freedom a topic of immense interest 
during the past decades. u > 

In nuclear physics, one of the typical examples of tunneling phenomena is heavy- 
ion fusion reaction at energies near and below the Coulomb barrier. 12 )' 13 ) Fusion 
is defined as a reaction in which two separate nuclei combine together to form a 
compound nucleus. In order for fusion reaction to take place, the relative motion be- 
tween the colliding nuclei has to overcome the Coulomb barrier formed by a strong 
cancellation between the long-ranged repulsive Coulomb and the short-ranged at- 
tractive nuclear forces (as a typical example, Fig.l shows the internucleus potential 
between 16 and 144 Sm nuclei as a function of the relative distance). Unless under 
extreme conditions, it is reasonable to assume that atomic nuclei are isolated sys- 
tems and the couplings to external environments can be neglected. Nevertheless, one 
can still consider intrinsic environments. The whole spectra of excited states of the 
target and projectile nuclei (as well as several sorts of nucleon transfer processes) 
are populated in a complex way during fusion reactions. They act as environments 
to which the relative motion between the colliding nuclei couples. In fact, it has 
by now been well established that cross sections of heavy-ion fusion reactions are 
substantially enhanced due to couplings to nuclear intrinsic degrees of freedom at 
energies below the Coulomb barrier as compared to the predictions of a simple po- 
tential model. 12 )~ 18 ) Heavy-ion subbarrier fusion reactions thus make good examples 
of environment- assisted tunneling phenomena. 

Theoretically the standard way to address the effects of the couplings between 
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the relative motion and nuclear intrinsic degrees of freedom on fusion reactions is 
to numerically solve the coupled-channels equations which include all the relevant 
channels. In the eigen-channel representation of coupled-channels equations, the 
channel coupling effects can be interpreted in terms of a distribution of fusion bar- 
riers. 13 )> 19 )~ 21 ) In this representation, the fusion cross section is given by a weighted 
sum of the fusion cross sections for each eigen-barrier. Those eigen-barriers lower 
than the original barrier are responsible for the enhancement of the fusion cross sec- 
tion at energies below the Coulomb barrier. Based on this idea, Rowley, Satchler, 
and Stelson have proposed a method to extract barrier distributions directly from 
experimental fusion excitation functions by taking the second derivative of the prod- 
uct of the fusion cross section and the center of mass energy Ea{ us with respect to E, 
i.e., d 2 (Ea ius )/dE 2 , 22 ) This method was tested against high precision experimental 
data of fusion cross sections soon after the method was proposed. 23 ) The extracted 
fusion barrier distributions were sensitive to the effects of channel-couplings and 
provided a much more apparent way of understanding their effects on the fusion 
process than the fusion excitation functions themselves. It is now well recognised 
that the barrier distribution approach is a standard tool for heavy-ion subbarrier 
fusion reactions. 13 ^' 18 ) 

The aim of this paper is to review theoretical aspects of heavy-ion subbarrier 
fusion reactions from the view point of quantum tunneling of composite particles. To 
this end, we mainly base our discussions on the coupled-channels approach. Earlier 
reviews on the subbarrier fusion reactions can be found in Refs. 12)-16). See also 
Refs. 24) and 25) for reviews on subbarrier fusion reactions of radioactive nuclei, 
and e.g., Refs. 26) and 27) for reviews on fusion reactions relevant to synthesis of 
superheavy elements, both of which we do not cover in this article. 

The paper is organized as follows. We will first discuss in the next section a po- 
tential model approach to heavy-ion fusion reactions. This is the simplest approach 
to fusion reaction, in which only elastic scattering and fusion are assumed to occur. 
This approach is adequate for light systems, but for fusion with a medium-heavy or 
heavy target nucleus the effects of nuclear excitations during fusion start playing an 
important role. In Sec. 3, we will discuss such nuclear structure effect on heavy- 
ion fusion reactions. To this end, we will introduce and detail the coupled-channels 
formalism which takes into account the inelastic scattering and transfer processes 
during fusion reactions. In Sec. 4, light will be shed on the fusion barrier distri- 
bution representation of fusion cross section defined as d 2 (Ea{ us )/dE 2 . It has been 
known that this approach is exact when the excitation energy of the intrinsic motion 
is zero, but we will demonstrate that one can generalize it unambiguously using the 
eigen-channel approach also to the case when the excitation energy is finite. In Sec. 
5, we will turn to a discussion on the present status of our understanding of deep 
subbarrier fusion reactions. At these energies, fusion cross sections have been shown 
to be suppressed compared to the values of the standard coupled-channels calcula- 
tions. This phenomenon may be related to dissipative quantum tunneling, that is, 
an irreversible coupling to intrinsic degrees of freedom. In Sec. 6, we discuss an ap- 
plication of the barrier distribution method to surface physics, more specifically, the 
effect of rotational excitations on a dissociative adsorption process of H2 molecules. 
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We then summarize the paper in Sec. 7. 

§2. One dimensional potential model 

2.1. Ion-ion potential 

Theoretically, the simplest approach to heavy-ion fusion reactions is to use the 
one dimensional potential model where both the projectile and the target are as- 
sumed to be structureless. A potential between the projectile and the target is given 
by a function of the relative distance r between them. It consists of two parts, that 
is, 

V(r) = V N (r)+V c (r), (2-1) 
where Vjv(r) is the nuclear potential, and Vc(r) is the Coulomb potential given by 

Voir) = (2-2) 
r 

in the outside region where the projectile and the target nuclei do not significantly 
overlap with each other. Figure 1 shows a typical potential V(r) for the s-wave 
scattering of the 16 + 144 Sm reaction. The dotted and the dashed lines are the 
nuclear and the Coulomb potentials, respectively, while the total potential V(r) is 
denoted by the solid line. One can see that a potential barrier appears due to a 
strong cancellation between the short-ranged attractive nuclear interaction and the 
long-ranged repulsive Coulomb force. This potential barrier is referred to as the 
Coulomb barrier and has to be overcome in order for the fusion reaction to take 
place. i?touch hi the figure is the touching radius, at which the projectile and the 
target nuclei begin overlapping considerably. One can see that the Coulomb barrier 
is located outside the touching radius. 

There are several ways to estimate the nuclear potential V/v(r). One standard 
method is to fold a nucleon-nucleon interaction with the projectile and the target 
densities. 28 ^ The direct part of the nuclear potential in this double folding procedure 
is given by 

Vn(t) = J dr 1 dr 2 v NN (r2-r 1 -r)pp(r 1 )p T (r 2 ), (2-3) 

where vnn is an effective nucleon-nucleon interaction, and pp and pr are the den- 
sities of the projectile and the target, respectively. The double-folding potential is 
in general a non-local potential due to the anti-symmetrization effect of nucleons. 
Usually, either a zero-range approximation 28 )' 29 ) or a local momentum approxima- 
tion 30 )~ 34 ) is employed in order to treat the non-locality of the potential. 

A phenomenological nuclear potential has also been employed. For instance, a 
Woods-Saxon form 

with 

V = 16vr 7 i?a, (2-5) 
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has been widely used, where the parameters were determined from a least-squares 
fit to the experimental data of heavy-ion elastic scattering. 35 )' 36 ) 

A nuclear potential so constructed has been successful in reproducing exper- 
imental angular distributions of elastic and inelastic scattering for many systems. 
Moreover, the empirical value of surface diffuseness parameter, a ~ 0.63 fm, is con- 
sistent with a double folding potential. Recently, a value of the surface diffuseness 
parameter has been determined unambiguously using heavy-ion quasi-elastic scatter- 
ing at deep subbarrier energies. 37 )' 38 ) It has been confirmed that the experimental 
data are consistent with a value around a ~ 0.63 fm. 38 )~ 41 ) 

In marked contrast, recent experimental data for heavy-ion subbarrier fusion re- 
actions suggest that a much larger value of diffuseness, ranging from 0.75 to 1.5 fm, 
is required to fit the data. 18 )' 42 ) -46 ) The Woods-Saxon potential which fits elastic 
scattering overestimates fusion cross sections at energies both above and below the 
Coulomb barrier, having an inconsistent energy dependence with the experimental 
fusion excitation function. A reason for the large discrepancies in diffuseness parame- 
ters extracted from scattering and fusion analyses has not yet been fully understood. 
However, it is probably the case that the double folding procedure is valid only in 
the surface region, while several dynamical effects come into play in the inner part 
where fusion is sensitive to. 

We summarize the relation between the surface diffuseness parameter a of a 
nuclear potential and the parameters of the Coulomb barrier, that is, the curvature, 
the barrier height, and the barrier position in Appendix A for an exponential and a 
Woods-Saxon potentials. 

2.2. Fusion cross sections 

In the potential model, the internucleus potential, V(r), is supplemented by an 
imaginary part, —iW(r), which mocks up the formation of a compound nucleus. 
One then solves the Schrodinger equation 



-2^ + yW 



iW(r) 



1(1 + l)h 2 
2fir 2 



E 



ui(r) = 0, 



(2-11) 



for each partial wave /, where /i is the reduced mass of the system, with the boundary 
conditions of 



ui(r) 
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Here, and ^ are the outgoing and the incoming Coulomb wave functions, 

respectively. S\ is the nuclear 5-matrix, and k = y / 2fiE/h 2 is the wave number 
associated with the energy E. 

If the imaginary part of the potential, W(r), is confined well inside the Coulomb 
barrier, one can regard the total absorption cross section as fusion cross section, i.e., 

<?UE) ~ a abs (£) = p ^(21 + 1) (1 - |^| 2 ) . (2-14) 

In heavy-ion fusion reactions, instead of imposing the regular boundary condition 
at the origin, Eq. (2-12), the so called incoming wave boundary condition (IWBC) is 
often applied without introducing the imaginary part of the potential, VF(r). 19 )' 47 ) 
Under the incoming wave boundary condition, the wave function has a form 



= Y k^r) ^ 6XP k ^ r ') dr '^j r < r abs, (2-15) 

at the distance smaller than the absorption radius r a b s , which is taken to be inside 
the Coulomb barrier. Here, ki{r) is the local wave number for the Z-th partial wave 
defined by 

«r,-^(.-V M -i2±^). 

The incoming wave boundary condition corresponds to the case where there is a 
strong absorption in the inner region so that the incoming flux never returns back. 
For heavy-ion fusion reactions, the final result is not sensitive to the choice of the 
absorption radius r a b s , and it is often taken to be at the pocket of the potential. 48 ) 
With the incoming wave boundary condition, 71 in Eq. (2-15) is interpreted as the 
transmission coefficient. Equation (2-14) is then transformed to 

i 

where Pi(E) is the penetrability for the Z-wave scattering defined as 

P l (E) = l-\S l \ 2 = \Ti\ 2 , (248) 

for the boundary conditions (2-13) and (2-15). The mean angular momentum of the 
compound nucleus is evaluated in a similar way as 

For a parabolic potential, Wong has derived an analytic expression for fusion cross 
sections, Eq. (2-17). 49 ) We will discuss it in Appendix B. 

The incoming wave boundary condition, Eq. (2-15), has two advantages over the 
regular boundary condition, Eq. (2-12). The first advantage is that the imaginary 
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Fig. 2. Comparison of experimental fusion cross sections for the 14 N+ 12 C system (the left panel) 
and 16 0+ 154 Sm system (the right panel) with results of the potential model calculations (the 
solid line). The height of the Coulomb barrier is around V b ~ 6.9 MeV and 59 MeV for 14 N+ 12 C 
and le O+ 154 Sm, respectively. The experimental data are taken from Refs. 50) and 18) for the 
14 N+ 12 C and the 16 0+ 154 Sm reactions, respectively. 

part of nuclear potential is not needed, and the number of adjustable parameters 
can be reduced. The second point is that the incoming wave boundary condition 
directly provides the penetrability Pi(E) = \Ti\ 2 and thus the round off error can 
be avoided in evaluating 1 — \Si\ 2 . This is a crucial point at energies well below the 
Coulomb barrier, where Si is close to unity. Notice that the incoming wave boundary 
condition does not necessarily correspond to the limit of W(r) — > oo, as the quantum 
reflection due to W{r) has to be neglected in order to realize it. The incoming wave 
boundary condition should thus be regarded as a different model from the regular 
boundary condition. 

2.3. Comparison with experimental data: success and failure of the potential model 

Let us now compare the one dimensional potential model for heavy-ion fusion 
reaction with experimental data. Figure 2 shows the experimental excitation func- 
tions of fusion cross section for 14 N+ 12 C (the left panel) and 16 0+ 154 Sm (the right 
panel) systems, as well as results of the potential model calculation (the solid lines). 
One can see that the potential model well reproduces the experimental data for the 
lighter system, 14 N + 12 C. On the other hand, the potential model apparently un- 
derestimates fusion cross sections for the heavier system, 16 + 154 Sm, although it 
reproduces the experimental data at energies above the Coulomb barrier, which is 
about 59 MeV for this system. 

In order to understand the origin for the failure of the potential model, Fig. 
3 shows the experimental fusion excitation functions for 16 + 144 ' 148 > 154 Sm reac- 
tions 18 ) and a comparison with the potential model (the solid line). To remove trivial 
target dependence, these are plotted as a function of center of mass energy relative 
to the barrier height for each system, and the fusion cross sections are divided by 
the geometrical factor, vri? 2 . With these prescriptions, the fusion cross sections for 
the different systems are matched with each other at energies above the Coulomb 
barrier, although one can also consider a more refined prescription. 51 )> 52 ) The barrier 
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Fig. 3. The experimental fusion cross sections for 16 0+ 144 ' 148 ' 154 Sm systems, taken from Ref. 18). 
In order to remove the trivial target dependence, the experimental fusion cross sections are 
divided by ivRf, where Rb is the position of the Coulomb barrier, and the energies are measured 
with respect to the barrier height, Vb, for each system. The solid line shows the result of the 
potential model calculation. 



height and the result of the potential model are obtained with the Akyiiz-Winther 
potential. 36 ^ One again observes that the experimental fusion cross sections are dras- 
tically enhanced at energies below the Coulomb barrier compared with the prediction 
of the potential model. Moreover, one also observes that the degree of enhancement 
of fusion cross section depends strongly on the target nucleus. That is, the enhance- 
ment for the 16 + 154 Sm system is order of magnitude, while that for the 16 + 
144 Sm system is about a factor of four at energies below the Coulomb barrier. This 
strong target dependence of fusion cross sections suggests that low-lying collective 
excitations play a role, as we will discuss in the next section. 

The inadequacy of the potential model has been demonstrated in a more trans- 
parent way by Balantekin et a/. 53 ) Within the semi-classical approximation, the 
penetrability for a one-dimensional barrier can be inverted to yield the barrier thick- 
ness. 54 ) Balantekin et al. applied such inversion formula directly to the experimental 
fusion cross sections in order to construct an effective internucleus potential. Assum- 
ing a one-dimensional energy-independent local potential, the resultant potentials 
were unphysically thin for heavy systems, often multi-valued potential. This result 
was confirmed also by the systematic study in Ref. 55) . These analyses have provided 
a clear evidence for the inadequacy of the one-dimensional barrier passing model for 
heavy- ion fusion reactions, and has triggered to develop the coupled-channels ap- 
proach, which we will discuss in the next section. 

In passing, we have recently applied the inversion procedure in a modified way to 
determine the lowest potential barrier among the distributed ones due to the effects 
of channel coupling. 56 ) The extracted potential for the 16 + 208 Pb scattering is 
well behaved, indicating that the channel coupling indeed plays an essential role in 
subbarrier fusion reactions. 
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Fig. 4. The experimental low-lying spectra of 144 > 148 > 154 Sm nuclei 



§3. Coupled-channels formalism for heavy-ion fusion reactions 

3.1. Nuclear structure effects on subbarrier fusion reactions 

The strong target dependence of subbarrier fusion cross sections shown in Fig. 3 
suggests that the enhancement of fusion cross sections is due to low-lying collective 
excitations of the colliding nuclei during fusion. The low-lying excited states in even- 
even nuclei are collective states, and strongly reflect the pairing correlation and shell 
structure. They have thus strongly coupled to the ground state, and also have a 
strong mass number and atomic number dependences. As an example, the low-lying 
spectra are shown in Fig. 4 for 144 > 148 > 154 Sm. The 144 Sm nucleus is close to the 
(sub-)shell closures (Z=64 and N = 82) and is characterized by a strong octupole 
vibration. 154 Sm, on the other hand, is a well deformed nucleus, and has a well 
developed ground state rotational band. 148 Sm is a transitional nucleus, and there 
exists a soft quadrupole vibration in the low-lying spectrum. One can clearly see 
that there is a strong correlation between the degree of enhancement of fusion cross 
sections shown in Fig. 3 and e.g., the energy of the first 2 + state. 

Besides the low-lying collective excitations, there are many other modes of exci- 
tations in atomic nuclei. Among them, non-collective excitations couple only weakly 
to the ground state and usually they do not affect in a significant way heavy-ion fu- 
sion reactions, even though the number of non-collective states is large. 57 ) Couplings 
to giant resonances are relatively strong due to their collective character. However, 
since their excitation energies are relatively large and also are smooth functions 
of mass number, 58 ) -60 ^ their effects can be effectively incorporated in a choice of 
internuclear potential through the adiabatic potential normalization (see the next 
section). 

The effect of rotational excitations of a heavy deformed nucleus can be easily 
taken into account using the orientation average formula. 17 )' 21 )' 49 )' 61 ) -63 ) For an ax- 
ially symmetric target nucleus, fusion cross sections are computed with this formula 
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Fig. 5. (The left panel) the orientation dependence of fusion potential for the 16 0+ 154 Sm reaction. 
The solid and the dashed lines are the potentials when the orientation of the deformed 154 Sm 
target is 6 — and 7r/2, respectively. The dotted line denotes the potential when the deformation 
of 154 Sm is not taken into account. (The right panel) the fusion cross sections for the 16 0+ 154 Sm 
reaction. The dashed line is the result of the potential model calculation shown in Fig. 3, while 
the solid line is obtained by taking into account the deformation of the 154 Sm nucleus with Eq. 
(3T). The experimental data are taken from Ref. 18). 



as, 



<Xfus(£) = / d(cos9)a {us (E;9), (3-1) 
J o 



where 9 is the angle between the symmetry axis and the beam direction. af us (E; 9) 
is a fusion cross section for a fixed orientation angle, 9. This is obtained with e.g., a 
deformed Woods-Saxon potential, 

Vo 

Fjv(r ' 9) = ~ 1 + exp[(r - R - Rt^Y^B) - R T ^Y m {e))/a\ ' (3 ' 2) 

which can be constructed by changing the target radius Rt in the Woods-Saxon 
potential, Eq. (24), to R T -»• R T (l + &*2o(0) + l3 A Y m (9)). See Ref. 64) for a 
recent application of this formula to fusion of massive systems, in which the formula 
is combined with classical Langevin calculations. 

The left panel of Fig. 5 shows the potential for the 16 0+ 154 Sm reaction obtained 
with the deformation parameters of /32=0.306 and j3 A = 0.05. The deformation of 
the Coulomb potential is also taken into account (see Sec. 3.4 for details). The 
solid line shows the potential for 9 = 0. For this orientation angle, the potential is 
lowered by the deformation effect as compared to the spherical potential shown by 
the dotted line, because the attractive nuclear interaction is active from relatively 
large values of r. The opposite happens when 9 = tt/2 as shown by the dashed 
line. The potential is distributed between the solid and the dashed lines accord- 
ing to the value of orientation angle, 9. The solid line in the right panel of Fig. 
5 shows the fusion cross sections obtained by averaging the contribution of all the 
orientation angles through Eq. (3-1). Since the tunneling probability has an ex- 
ponentially strong dependence on the barrier height, the fusion cross sections are 
significantly enhanced for those orientations which yield a lower barrier than the 
spherical case. It is remarkable that this simple calculation accounts well for the 
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experimental enhancement of fusion cross sections at subbarrier energies. Evidently, 
nuclear structure effects significantly enhance fusion cross sections at energies below 
the Coulomb barrier, which make fusion reactions an interesting probe for nuclear 
structure. 

3.2. Coupled- channels equations with full angular momentum coupling 

The nuclear structure effects can be taken into account in a more quantal way 
using the coupled-channels method. In order to formulate the coupled-channels 
method, consider a collision between two nuclei in the presence of the coupling of 
the relative motion, r = (r, f), to a nuclear intrinsic motion £. We assume the 
following Hamiltonian for this system, 

H(r, = -|V + V(r) + H (0 + V coup (r, £), (3-3) 

where i?o(£) and V coup (r,^) are the intrinsic and the coupling Hamiltonians, respec- 
tively. In general the intrinsic degree of freedom £ has a finite spin. We therefore 
expand the coupling Hamiltonian in multipoles as 

Wr,£) = ^fx(r)Y x (r) -T A (£). (3-4) 

A>0 

Here Y\(r) are the spherical harmonics and are spherical tensors constructed 

from the intrinsic coordinate. The dot indicates a scalar product. The sum is taken 
over all values of A except for A = 0, which is already included in the bare potential, 
V(r). 

For a given total angular momentum J and its z component M, one can define 
the channel wave functions as 

(H\(alI)JM) = ^ (lm l Im I \JM)Y lmi (r) VaImi ^), (3-5) 

mi,mi 

where I and / are the orbital and the intrinsic angular momenta, respectively. 
(£) are the wave functions of the intrinsic motion which obey 

(0- (3-6) 

Here, a denotes any quantum number besides the angular momentum. Expanding 
the total wave function with the channel wave functions as 

*j(r,0 = £ ^l(rC\(alI)JM), (3-7) 

ce,l,I 

the coupled-channels equations for u J alI {r) read 

J aii(r)+ Yl V ali;<x'l>i>(r)u J a , llII (r) = 0, 



h 2 d 2 i(l + i)h 2 , , ^ 



(3-8) 
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where the coupling matrix elements Vaira'i'i'( r ) are gi ven as 65 -* 

Vama'i'Ar) = ((alI)JM\V coup (r,0\(al'l')JM), (3-9) 



= E(-) / " 7 ' +z ' +J /AW(/nniiO(^iiTAiia7') 

A 

V(2/ + l)(2/ + l)| f J ^ }. (3-10) 



A 

X- 



Notice that these matrix elements are independent of M. 

For the sake of simplicity of the notation, in the following let us introduce a 
simplified notation, n = {a, I, I}, and suppress the index J. The coupled-channels 
equation (3-8) then reads, 



Mr) = 0. (3-11) 



These coupled-channels equations are solved with the incoming wave boundary con- 
ditions of 



Mr) ~ J t¥t 7~rl exp ( -i I k n (r')dr' ) r < r abs , (3-12) 



M r ) 



r abs 



= H ( ~\k n r)S n , ni - J^si H H^\k n r) r oo, (3-13) 

V 

where rij denotes the entrance channel. The local wave number fc n (r) is defined by 



Mr) ^|( £ _ e „_«k±^_ FW ), ,3,4, 

whereas fc„ = fc ra (r = oo) = y / 2/j>(E — e n )/h 2 . Once the transmission coefficients 
are obtained, the inclusive penetrability of the Coulomb potential barrier is 
given by 

Pj(E) = ^\T n J n f. (3-15) 

n 

The fusion cross section is then given by 

<tUE) = ^ J>J + l)Pj(^), (3-16) 

where we have assumed that the initial intrinsic state has spin zero, Ii = 0. This 
equation for fusion cross section is similar to Eq. (2-17) except that the penetrability 
Pj{E) is now influenced by the channel coupling effects. 

3.3. I so- centrifugal approximation 

The full coupled-channels calculations (3-11) quickly become intricate if many 
physical channels are included. The dimension of the resulting coupled-channels 
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problem is in general too large for practical purposes. For this reason, the iso- 
centrifugal approximation, which is sometimes referred to as the no-Coriolis approxi- 
mation or the rotating frame approximation, has often been introduced. 21 )' 48 )' 63 )' 66 ) -70 ) 
In the iso-centrifugal approximation to the coupled-channels equations, Eq. (3-11), 
one first replaces the angular momentum of the relative motion in each channel by 
the total angular momentum J, that is, 

i n {i n + \)h 2 ^ J(J + i)ft 2 

~^ 2 (3 ' 17) 

This corresponds to assuming that the change in the orbital angular momentum due 
to the excitation of the intrinsic degree of freedom is negligible. Introducing the 
weighted average wave function 

Mr) = (-Y 2(J0/0|Z0)«,j(r), (3-18) 
l 

where we have suppressed the index a for simplicity, and using the relation 

5Z(-) l ' + " 7+A V^r=FT | J x \, \, }(fflA0|/'0)(J0/0|ffl) = J^^ JOJ'OlZ'OK/'OAOl/O), 

(3-19) 

one finds that the wave function uj{r) obeys the reduced coupled-channels equations, 

/ h 2 d 2 J(J + i)h 2 , . „ \ , , v-v- /2A + 1 , ,. Im . . , , 

(3-20) 

These are nothing but the coupled-channels equations for a spin-zero system with 
the interaction Hamiltonian given by 



^coup = X>(r)n(r = 0) -T X = Y J ]J^^fx(r)T X o. (3-21) 

A A 

In solving the reduced coupled-channels equations, similar boundary conditions are 
imposed for uj as those for uu, 



" /(r) ~Vfc|)^ eXP (" ! i Mr')dr^j r<r abs , (3-22) 

= H^ik^Su, - ^S-ljM^ikjr) r -> oo, (3-23) 

where ki and kj(r) are defined in the same way as in Eq. (3T4). The fusion cross 
section is then given by Eq. (3-16) with the penetrability of 

P J (E) = J2\T/ Ii \ 2 . (3-24) 
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Since the reduced coupled-channels equations in the iso-centrifugal approxima- 
tion are equivalent to the coupled-channels equations with a spin-zero intrinsic mo- 
tion, the complicated angular momentum couplings disappear. A remarkable fact is 
that the dimension of the coupled-channels equations is drastically reduced in this 
approximation. For example, if one includes four intrinsic states with 2 + , 4 + , 6 + , 
and 8 + together with the ground state in the coupled-channels equations, the original 
coupled-channels have 25 dimensions for J > 8, while the dimension is reduced to 5 
in the iso-centrifugal approximation. The validity of the iso-centrifugal approxima- 
tion has been well tested for heavy- ion fusion reactions, and it has been concluded 
that the iso-centrifugal approximation leads to negligible errors in calculating fusion 
cross sections. 63 )> 67 ) 

3.4. Coupling to low-lying collective states 

3.4.1. Vibrational coupling 

Let us now discuss the explicit form of the coupling Hamiltonian V coup for heavy- 
ion fusion reactions. We first consider couplings of the relative motion to the 2 A - 
pole surface vibration of a target nucleus. In the geometrical model of Bohr and 
Mottelson, the radius of the vibrating target is parameterized as 

R(0, ^ = R T (^l + ^2 CLxJT^ip, ^ , (3-25) 

where Rt is the equivalent sharp surface radius and is the surface coordinate of 
the target nucleus. To the lowest order, the surface oscillation is approximated by a 
harmonic oscillator and the Hamiltonian for the intrinsic motion is given by 



o = tiw x (j2 + ^f^J ■ (3-26) 

Here hui\ is the oscillator quanta and and are the phonon creation and 
annihilation operators, respectively. The surface coordinate is related to the 
phonon creation and annihilation operators by 

= a « + (-fax,) = -j^= (4 M + (-Tax,) , (3-27) 

where «o = /?a/v / 2A~+T is the amplitude of the zero point motion. 58 ** The defor- 
mation parameter fix can be estimated from the experimental transition probability 
using (see Eq. (3-34) below) 



ftv = * J*£»l. (3-28) 
3Z T R*\ e 2 



The surface vibration of the target nucleus modifies both the nuclear and the 
Coulomb interactions between the colliding nuclei. In the collective model, the nu- 
clear interaction is assumed to be a function of the separation distance between the 
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vibrating surfaces of the colliding nuclei, and thus is given as 

V {N \r,a^) = V N [r- R T £ a^Y^r) j . (3-29) 

If the amplitude of the zero point motion of the vibration is small, one can expand 
this equation in terms of and keep only the linear term, 

yW(r,a v ) = V N (r)-R T ^1^2a x ,Y^(r). (3-30) 

This approximation is called the linear coupling approximation. The first term of 
the right hand side (r.h.s.) of Eq. (3-30) is the bare nuclear potential in the absence 
of the coupling, while the second term is the nuclear component of the coupling 
Hamiltonian. Even though the linear coupling approximation does not work well for 
heavy-ion fusion reactions, 48 - 1 ' 71 -' we employ it in this subsection in order to illustrate 
the coupling scheme. In Sec. 3.5, we will discuss how the higher order terms can be 
taken into account in the coupling matrix. 

The Coulomb component of the coupling Hamiltonian is evaluated as follows. 
The Coulomb potential between the spherical projectile and the vibrating target is 
given by 

\r / \ f, ,Z P Z T e 2 , Z P Z T e 2 ^ ^ 4-irZ P e 1 

(3-31) 

where pr is the charge density of the target nucleus and Qav the electric multipole 
operator defined by 

Qav = / drZ T ep T (r)r x 'Y yil ,(r). (3-32) 

The first term of the r.h.s. of Eq. (3-31) is the bare Coulomb interaction, and the 
second term is the Coulomb component of the coupling Hamiltonian. In obtaining 
Eq. (3-31), we have used the formula 



4"7T 



„A' 



' I A'/i' ^ r > 

and have assumed that the relative coordinate r is larger than the charge radius of 
the target nucleus. If we assume a sharp matter distribution for the target nucleus, 
the electric multipole operator is given by 

3e 

Q*v' = -^ z TRr a \^\ tJ ,,\' l x', (3-34) 

up to the first order in the surface coordinate a\^. 

By combining Eqs. (3-30), (3-31), and (3-34), the coupling Hamiltonian is ex- 
pressed by 

V coup (r, a x ) = fx(r) ^ a Xfl Y^(r), (3-35) 
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up to the first order of ax^. Here, f\(r) is the coupling form factor, given by 

/A(r) = _^ + _A_ ZpZre ^, (3.36) 

where the first and the second terms are the nuclear and the Coulomb coupling form 
factors, respectively. Transforming to the rotating frame, the coupling Hamiltonian 
used in the iso-centrifugal approximation is then given by (see Eq. (3-21)), 



Vcoupir, a xo ) = \j 2X ^ 1 f\{r)a xo = -j=f\(r) (a{ + a A0 ) . (3-37) 

Notice that the coupling form factor f\ has the value 

, Z P Z T e* ( 3 R T \ ( 

fx{Rb) ~ if ~ bTJ ■ (3 ' 38) 

at the position of the bare Coulomb barrier, R^, and the coupling strength is ap- 
proximately proportional to the charge product of the colliding nuclei. 

In the previous subsection, we showed that the iso-centrifugal approximation 
drastically reduces the dimension of the coupled-channels equations. A further re- 
duction can be achieved by introducing effective multi-phonon channels. 66 )' 69 - ) In 
general, the multi-phonon states of the vibrator have several levels, which are dis- 
tinguished from each other by the angular momentum and the seniority. 58 -* For ex- 
ample, for the quadrupole surface vibrations, the two-phonon state has three levels 
(0 + , 2 + , 4 + ), which are degenerate in energy in the harmonic limit. The one-phonon 
state, \2~l ) = a^o 1 0), couples only to a particular combination of these triplet states, 

|2>= Y, <2020|/0>|/0> = ^(4 ) 2 |0>. (3-39) 

7=0,2,4 

It is thus sufficient to include this single state in the calculations, instead of three 
triplet states. In the same way, one can introduce the n-phonon channel for a 
multipolarity A as 

See Appendix C for the case of two different vibrational modes of excitation {e.g., a 
quadrupole and an octupole vibrations). 

If one truncates the phonon space up to the two-phonon state, the corresponding 
coupling matrix is then given by 



^coup = F(r) Hw x V2F(r) , (3-41) 




where F(r) is defined as fi\f\{r)/y/ ^tt. 

The effects of deviations from the harmonic oscillator limit presented in this 
subsection on subbarrier fusion reactions have been discussed in Refs. 72) and 73). 
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3.4.2. Rotational coupling 

We next consider couplings to the ground rotational band of a deformed target. 
To this end, it is convenient to transform to the body fixed frame where the z axis 
is along the orientation of the deformed target. The surface coordinate is then 
transformed to 

a\^ = ^2D* lfl ((j) d ,6 d ,Xd)a\ fl >, (3-42) 
uf 

where 4> d ,Qd, an d Xd are the Euler angles which specify the body-fixed frame, thus 
the orientation of the target. If we are particularly interested in the quadrupolc 
deformation (A=2), the surface coordinates in the body fixed frame are expressed as 

a 2 o = #2 cos 7, (3-43) 

022 = a 2 _ 2 = ^/?2sin7, (3-44) 
v2 

021 = a 2 _i = 0. (3-45) 

If we further assume that the deformation is axial symmetric (i.e., 7 = 0), the 
coupling Hamiltonian for the rotational coupling reads (see Eq. (3-35)) 

V coup (r,9 d ,<t> d ) = f2{r)Y^^\j^ Y ^ 9 dAd)Y^{r). (3-46) 
In order to obtain this equation, we have used the relation 



DU(^ 0,x) = \j ^Ty Y *lm(0, 0). (3-47) 

The coupling Hamiltonian in the rotating frame is thus given by 

V coup (r,6) = / 2 (r)/3 2 y 2O (0), (3-48) 

where 9 is the angle between (9 d , (fi d ) and f , that is, the direction of the orientation of 
the target measured from the direction of the relative motion between the colliding 
nuclei. Since the wave function for the |I0) state in the ground rotational band is 
given by |/0) = \Y[q), the corresponding coupling matrix is given by 



Koup = f2(r)/3 2 (Y ro \Y 20 \Y IO ) = 



( F(r) \ 

F(r) e 2 + ^F(r) f F(r) 

V f F(r) fe 2 + 2 -^F(r) j 



(3-49) 

when the rotational band is truncated at the first 4 + state. Here, e 2 is the excitation 
energy of the first 2 + state, and F(r) is defined as ^.MO/V^vr as in Eq. (3-41). 

One of the main differences between the vibrational (3-41) and the rotational 
(3-49) couplings is that the latter has a diagonal component which is proportional to 
the deformation parameter /3 2 . The diagonal component in the rotational coupling 
is referred to as the reorientation effect and has been used in the Coulomb excitation 
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technique to determine the sign of the deformation parameter. 74 ) Notice that the 
results of the coupled-channels calculations are independent of the sign of j3 2 for the 
vibrational coupling. 

The effects of the 7 deformation on subbarrier fusion were studied in Ref. 75). 
If there is a finite 7 deformation, the coupling Hamiltonian in the rotating frame 
becomes 

Kou P (r, 9, 0) = f 2 (r) (fa cos 7 y 20 (#) + -j=/3 2 sin 7 (Y 22 (9, 0) + Y 2 „ 2 (9, 0)) 

(3-50) 

Higher order deformations can also be taken into account in a similar way as the 
quadrupole deformation. For example, if there is an axial symmetric hexadecapole 
deformation in addition to quadrupole deformation, the coupling Hamiltonian reads 

V COVLp (r, 9) = f 2 (r)foY 2Q {6) + f 4 (r)p 4 Y 40 (9), (3-51) 

where (3^ is the hexadecapole deformation parameter. 

3.5. All order couplings 

In the previous subsection, for simplicity, we have used the linear coupling ap- 
proximation and expanded the coupling Hamiltonian in terms of the deformation 
parameter. However, it has been shown that the higher order terms play an impor- 
tant role in heavy-ion subbarrier fusion reactions. 48 **' 71 ) ,76 )~ 79 ) These higher order 
terms can be evaluated as follows. 48 ) If we employ the Woods-Saxon potential, Eq. 
(2-4), the nuclear coupling Hamiltonian can be generated by changing the target 
radius in the potential to a dynamical operator 

R ^R + O, (3-52) 

that is, 

V N (r) -> V N (r, 6) = -— V ° . ■ (3-53) 

1 + exp((r - R - 0)/a) 

For the vibrational coupling, the operator O is given by (see Eq. (3-37)), 

6 = ^=R T (a[ + a X0 ), (3-54) 



while for the rotational coupling it is given by (see Eqs. (3-2) and and (3-48)), 

6 = (3 2 R T Y 20 (9) + p A R T Y m {6). (3-55) 

The matrix elements of the coupling Hamiltonian can be easily obtained using a 
matrix algebra. 80 ) In this algebra, one first looks for the eigenvalues and eigenvectors 
of the operator O which satisfies 

6\a) = \ a \a). (3-56) 

This is done by numerically diagonalising the matrix O, whose elements are given 
by 

Onm = —7=RT{VrnS njm -i + \/™Wn+l)- ( 3 '57) 
V47T 
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Fig. 6. The barrier penetrability for a two-level problem as a function of energy E in the linear 
(the left panel) and the logarithmic (the right panel) scales. The solid and the dashed lines are 
the exact solution and the WKB approximation, respectively. 



for the vibrational case, and 



/ 5(2/ + l)(2I' + l) (I 2 V ^ 2 
Ojp-\I ^ foR T [ Q Q Q 



/ 9(2j+i)(2/' + iy fl (i 4 r y 

+ V 4^ P* R t[ o o ) ■ (3-58) 

for the rotational case. The nuclear coupling matrix elements are then evaluated as 

V<t% ] = (n\V N (r,6)\m) - V N (r)S n , m , 

= y^(n\a)(a\m)V N (r, X a ) - V N (r)5 n)m . (3-59) 

a 

The last term in this equation is included to avoid the double counting of the diagonal 
component. 

The computer code CCFULL has been written with this scheme, 48 ) and has been 
used in analysing recent experimental fusion cross sections for many systems. CCFULL 
also includes the second order terms in the Coulomb coupling for the rotational case, 
while it uses the linear coupling approximation for the Coulomb coupling in the 
vibrational case. 48 ) 

3.6. WKB approximation for multi-channel penetrability 

Whereas the coupled-channels equations, Eq. (3-20), can be numerically solved 
e.g., with the computer code CCFULL once the coupling Hamiltonian has been set 
up, it is always useful to have an approximate solution. In the next section, we will 
discuss the limit of zero excitation energy for intrinsic degrees of freedom, in which 
the coupled-channels equations are decoupled. In this subsection, on the other hand, 
we discuss another approximate solution based on the semiclassical approximation. 

The penetrability in the WKB approximation is well known for a one dimensional 
potential V(x) and is given by, 



P(E) = exp 



2£dx'^(V(x')-E) 



(3-60) 
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where xo and x\ are the inner and the outer turning points satisfying V(xq) = 
V(xi) = E, respectively. One can also introduce the uniform approximation to take 
into account the multiple reflection under the barrier, and obtain a formula which is 
valid at all energies from below to above the barrier, 81 )~ 85 ) 

P(E) = = 1 = -. (3-61) 



1 + exp 



(V(x') - E) 



It has been shown in Ref. 86) that one can generalize the primitive WKB formula 
(3-60) to a multi-channel problem as, 



n 



,iq(xi)Ax 



(3-62) 



where q(x) = [2/x(£- W(x))/ft 2 ] 1/2 with W nm {x) = (n\V(x)+H (0+V coup (x,0\m) 
(see Eq. (3-3)). Here we have discretized the coordinate x with a mesh spacing of 
Ax. For a single channel problem, Eq. (3-62) is reduced to Eq. (3-60). 

Figure 6 shows the result of the multi-channel WKB approximation for a two- 
level problem given by 



W(x) 



with 



V(x) 
F(x) 



F(x) 
V(x) + e 



V(x) 



+ F(x) 



+ 



V(x) = V e- x2 / 2s ' 2 , F(x) = F e- x2/2s f 





e 

(3-63) 
(3-64) 



The parameters are chosen following Ref. 19) to be Vo=100 MeV, Fo=3 MeV, and 
s = Sf =3 fm, which mimic the fusion reaction between two 58 Ni nuclei. The 
excitation energy e and the mass /x are taken to be 2 MeV and 29rriN, respectively, 
where tun is the nucleon mass. It is remarkable that the WKB formula (3-62) 
reproduces almost perfectly the exact solution at energies well below the barrier. 
The WKB formula breaks down at energies around the barrier, as in the single- 
channel problem. 

The figure also suggests that the penetrability is given by a weighted sum of two 
penetrabilities, 

P(E) = w 1 P{E-\ 1 {x)) + w 2 P{E-\ 2 {x)), (3-65) 

where A«(x) are the eigen-potentials, Xi(x) = V(x) + [e± ^ e 2 + 4F(x) 2 ]/2, obtained 
by diagonalizing the matrix W(x), (3-63). We will discuss this point in the next 
section. 



§4. Barrier distribution representation of mult i- channel penetrability 

4.1. Sudden tunneling limit and barrier distribution 

In the limit of vanishing excitation energy for the intrinsic motion (i.e., in the 
limit of ej — > 0), the reduced coupled-channels equations (3-20) are completely de- 
coupled. This limit corresponds to the case where the tunneling occurs much faster 
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than the intrinsic motion, and thus is referred to as the sudden tunneling limit. In 
this limit, the coupling matrix defined as 



Vir = + y * * 1 f\{r){yio\Txo\vm) (4-1) 

can be diagonalized independently of r (for simplicity we consider only a single value 
of A). See also Eq. (3-63). It is then easy to prove that the fusion cross section is 
given as a weighted sum of the cross sections for uncoupled eigenchannels, 21 )' 87 ) 

aUE) = ^ w -4:l(E), (4-2) 

a 

where o'^(E) is the fusion cross section for a potential in the eigenchannel a, i.e., 
Va(r) = V(r) + A a (r). The same relation holds also for the quasi-elastic scatter- 
ing^), 87), 88) jj ere; A a (r) is the eigenvalue of the coupling matrix (4-1) (when e/ 
is zero, A Q (r) is simply given by A a • f\(r)). The weight factor w a is given by 
Wa = | t/oci | 2 , where U is the unitary matrix which diagonalizes Eq. (4-1). Note that 
the unitarity of the matrix U leads to the relation that the sum of all the weight 
factors, ^2 a w a , is unity. 21 ) 

The resultant formula (4-2) in the sudden tunneling limit can be interpreted in 
the following way. In the absence of the coupling, the incident particle encounters 
only the single potential barrier, V(r). When the coupling is turned on, the bare 
potential splits into many barriers. Some of them are lower than the bare potential 
and some of them higher. In this picture, the potential barriers are distributed with 
appropriate weight factors, w a . 

The orientation average formula discussed in Sec. 3.1 (see Eq. (3-1)) for a 
deformed target nucleus can also be obtained from the coupled-channels equations 
by taking the sudden tunneling limit. 21 ) To show this, first notice that the coupling 
Hamiltonian is diagonal with respect to the orientation angle, 9. If all the members 
of the rotational band are included in the coupled-channels equations, the eigenstates 
of the coupling Hamiltonian matrix then become the same as the angle vector \8) 
with the eigenvalue given by the deformed Woods-Saxon potential, Eq.(3-2). 21 )' 89 )' 90 ) 
The weight factor in this case is simply given by w{6) = |(#|^/=o)| 2 = |loo(#)| 2 - 

The physical interpretation of the orientation average formula is that fusion 
reaction takes place so suddenly that the orientation angle is fixed during the fusion 
reaction. This is justified because the first 2 + state of a heavy deformed nucleus 
is small (see Fig. 4), corresponding to a large moment of inertia for the rotational 
motion. As the orientation angles are distributed according to the wave function 
for the ground state, the fusion cross section can be computed first by fixing the 
orientation angle and then averaging over the orientation angle with the appropriate 
weight factor, w{9). The applicability of this formula has been investigated in Ref. 
62) in the reactions of 154 Sm target with various projectiles ranging from 12 C to 
40 Ar. It has been shown that the formula works well, although the agreement with 
the exact coupled-channels calculations which take into account the finite excitation 
energy of the rotational excitation becomes slightly worse for a large value of charge 
product of the projectile and the target nuclei. 
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Fig. 7. The product of energy E and the fusion cross section a{ us , Eaf UB , for the 16 0+ 144 Sm reaction 
obtained with the potential model (the left panel). The middle and the right panels show the 
first and the second energy derivatives of Eaf us , respectively. 



4.2. Fusion barrier distribution 

Rowley, Satchler, and Stelson have proposed a method to extract, directly from 
the experimental fusion cross sections, a way how the barriers are distributed. 13 )> 22 ) 
In order to illustrate the method, let us first discuss the classical fusion cross section 
given by, 

at(E) = ttRI (l - ^) 6(E - V b ). (4-3) 

From this expression, it is clear that the first derivative of Eaf^ s is proportional to 
the classical penetrability for a 1-dimensional barrier of height V b , 

-^[Eat(E)} = nR 2 0{E - V b ) = nR 2 P d (E), (4-4) 
and the second derivative to a delta function, 

"' 2 [Eo-jl{E)}=vRl5{E-V b ). (4-5) 



dE 2 



In quantum mechanics, the tunneling effect smears the delta function in Eq. 
(4-5). As we have noted in Sec. 2.2, an analytic formula for the fusion cross section 
can be obtained if one approximates the Coulomb barrier by an inverse parabola, 
see Eq. (B-5) in Appendix B. Again, the first derivative of Eo~f us is proportional to 
the s-wave penetrability for a parabolic barrier, 

±{Eo ins {E)\ = ttR 2 — I = irRl P(E), (4-6) 

dE 1 + exp [-%{E- V b )\ 

and the second derivative is proportional to the derivative of the s-wave penetrability, 

d 2 /^i 2vr e x „?dP(E) , . 

As shown in Fig. 7, this function has the following properties: i) it is symmetric 
around E = V b , ii) it is centered at E = V b , iii) its integral over E is irR b , and iv) it 
has a relatively narrow width of around ln(3 + y8)hf2/Tr ~ 0.56^1?. 
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Fig. 8. (The left panel) The fusion barrier distribution D lus (E) = d 2 (Etj{ ua )/dE 2 for the 
16 0+ 154 Sm reaction. 18 ' The solid line is obtained with the orientation average formula, which 
corresponds to the solid line in Fig. 5. The dashed lines indicate the contributions from six 
individual eigenbarriers (i.e., orientation angles). (The right panel) The fusion barrier distri- 
bution for the le O+ 144 Sm reaction. 18 ' The solid line shows the result of the coupled-channels 

144 72^ 73^ 

calculations which take into account the anharmonic double phonon excitations of Sm. ' 
In the presence of channel couplings, Eq. (4-2) immediately leads to 

Aus = ^[E*am = j>« ^[E^m- (4-8) 

a 

This function has been referred to as fusion barrier distribution. As an example, 
the left panel of Fig. 8 shows the barrier distribution for the 16 0+ 154 Sm reaction, 
whose fusion cross sections have been already shown in Fig. 5. We replace the 
integral in Eq. (3-1) with the (i max + 2)-point Gauss quadrature with J max =10. This 
corresponds to taking 6 different orientation angles. 21 ) The contributions from each 
eigenbarrier are shown by the dashed line in Fig. 8. The solid line is the sum of all 
the contributions, which is compared with the experimental data. 18 ) One can see 
that the calculation well reproduces the experimental data. Moreover, this analysis 
suggests that 154 Sm is a prolately deformed nucleus, since if it were an oblate nucleus, 
then lower potential barriers would have larger weights, and -D ms would be larger for 
smaller E, in contradiction to the experimentally observed barrier distribution. 13 ) 

The fusion barrier distribution has been extracted for many systems, see Ref. 
13) and references therein. The extracted barrier distributions were shown to be 
sensitive to the effects of channel-couplings and have provided a much more apparent 
way of understanding their effects on the fusion process than the fusion excitation 
functions themselves. These experimental data have thus enabled a detailed study 
of the effects of nuclear intrinsic excitations on fusion reactions, and have generated 
a renewed interest in heavy-ion subbarrier fusion reactions. An important point is 
that the nature of sub-barrier fusion reactions as a tunneling process exponentially 
amplifies the effects of detailed nuclear structure. Fusion barrier distribution makes 
this effect even more visible by plotting in the linear scale. The sub-barrier fusion 
reactions thus offer a novel way of nuclear spectroscopy which could be called a 
tunneling assisted nuclear spectroscopy. As an example, we mention that it was 
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recently applied to elucidate the shape transition and the shape coexistence of Ge 
isotopes. 91 ^' 92 ) It is worthwhile to mention also that the method of the barrier 
distribution has been successfully applied to heavy- ion quasi-elastic scattering. 63 )' 93 ) 

4.3. Eigenchannel representation 

As we have discussed in the previous subsection, the barrier distribution repre- 
sentation, that is, the second derivative of E<jf us , has a clear physical meaning only 
if the excitation energy of the intrinsic motion is zero. The concept holds only ap- 
proximately when the excitation energy is finite. Nonetheless, this analysis has been 
successfully applied to systems with relatively large excitation energies. 18 )' 79 )' 94 ) For 
example, the second derivative of Eaf us for 16 + 144 Sm fusion reaction has a clear 
double-peaked structure (see the right panel of Fig. 8). 18 )' 94 ) The coupled-channels 
calculation also yields a similar double-peaked structure of the fusion barrier distri- 
bution, and this structure has been interpreted in terms of the anharmonic octupole 
phonon excitations in 144 Sm, 72 )' 73 ) whose excitation energy is 1.8 MeV for the first 
3~ state. Also the analysis of the fusion reaction between 58 Ni and 60 Ni, where the 
excitation energies of quadrupole phonon states are 1.45 and 1.33 MeV, respectively, 
shows that the barrier distribution representation depends strongly on the number of 
phonon states included in coupled-channels calculations. 79 ) These analyses suggest 
that the representation of fusion process in terms of the second derivative of Eaf us is 
a powerful method to study the details of the effects of nuclear structure, irrespective 
of the excitation energy of the intrinsic motion. 

When the excitation energy of the intrinsic motion is finite, the barrier distribu- 
tion is defined in terms of the eigen-channels. To illustrate it, first notice that Eq. 
(3-15) can be expressed as 

P(E) = (TtT) nini , (4-9) 

using the completeness of the channels n (we have suppressed the index J) . We then 
introduce the eigenfunctions of the Hermitian operator T^T as, 

(TtT)|^) = 7fc |&>. (4-10) 

Using this basis, the penetrability is given by 

P{E) = Y,m\n t )\ 2 - lk . (4-11) 

k 

When the excitation energies e n are all zero, as we have discussed in Sec. 4.1, one 
can diagonalize the coupling matrix V nn > (r) with the basis set which is independent 
of the radial coordinate r. In this case, the matrix T is diagonal on this basis, and 
the weight factor \(cj)k\ni}\ 2 is independent of E. Eq. (4-11) is a generalization of 
this scheme, which is applicable also when the excitation energies are non-zero. 

Figure 9 shows the two eigenvalues 7^ and the corresponding weight factors 
\{ ( t ) k\ n i)\ 2 as a function of E for a single-phonon coupling calculation for the s-wave 
16 0+ 144 Sm reaction. To this end, we have taken into account couplings to the single 
octupole phonon state in 144 Sm at 1.81 MeV with the deformation parameter of ^3 
= 0.205. The total probability P(E), and the penetrability of the two eigenbarriers, 
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Fig. 9. (The left panel) The s-wave penetrabilities for the 16 0+ 144 Sm reaction. The dotted line 
is obtained with the coupled-channels calculations with a single octupole phonon excitation in 
144 Sm at 1.81 MeV with fls = 0.205. The solid lines show the eigenvalues of the square of the 
transmission matrix, T^T, defined by Eq. (4T0). The dashed lines denote the penetrabilities of 
the eigen barriers constructed by diagonalizing the coupling matrix at each r. (The right panel) 
The weight factors \{4>k\ni)\ 2 defined in Eq. (4-11) as a function of energy. 



obtained by diagonalizing the coupling matrix V nn >(r) at each r, are also shown in 
the left panel of the figure by the dotted and the dashed lines, respectively. One can 
see that the two eigenvalues 7^ approximately correspond to the penetrability of the 
eigenbarriers, and thus the factors | |nj)| 2 can be interpreted as the weight factors 
for each eigenbarrier. This implies that the fusion cross sections are still given by Eq. 
(4-2) even when the excitation energy is finite, except that the eigenbarriers are now 
constructed by diagonalizing the coupling matrix at each r. The weight factors do 
not vary strongly as a function of energy, suggesting that the concept of the fusion 
barrier distribution still holds to a good approximation even when the excitation 
energy of the intrinsic motion is finite. We have reached the same conclusion already 
in Ref. 95) using a different method from the one in this subsection. In contrast 
to the method in Ref. 95), the method in this subsection is more general since the 
applicability is not restricted to a two-level problem. 

4.4. Adiabatic potential renormalization 

Given that the concept of fusion barrier distribution still holds even with a finite 
excitation energy, it is interesting to investigate how the fusion barrier distribution 
evolves as the excitation energy is varied. To this end, we carry out the coupled- 
channels calculations for the 16 0+ 144 Sm reaction by taking into account the single 
octupole phonon excitation in 144 Sm. The solid line in Fig. 10 (a) shows the fusion 
barrier distribution Df us when the excitation energy of the octupole vibration, E s - , 
is set to be zero. For comparison, the figure also shows the result of no-coupling 
calculation by the dashed line. In this case, the original single barrier splits into two 
eigenbarriers with equal weight, one corresponds to the effective channel |0 + ) + |3~) 
and the other corresponds to |0 + ) — |3 _ ). The fusion barrier distribution is slightly 
asymmetric since the barrier positions, Rt>, are different between the two effective 
channels (see Eq. (4-7)). 

Figure 10(b) corresponds to the physical case of Eg- = 1.81 MeV. In this case, 
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Fig. 10. The fusion barrier distribution D{ us for the 16 0+ 144 Sm reaction with several values of 
excitation energies, E A - , of the octupole vibration in 144 Sm. The solid lines are the results of 
the coupled-channels calculations which take into account the single octupole phonon excitation 
in 144 Sm, while the dashed lines are obtained without taking into account the channel coupling 
effect. The curvature hQ of the Coulomb barrier is 4.25 MeV in these calculations. 




the barrier distribution still has a clear double peaked structure as in the experimen- 
tal data, 18 )' 94 ) but the lower energy barrier acquires more weight and the barrier 
distribution is highly asymmetric. The effective channels are now a|0 + ) + /3|3~) (the 
lower energy barrier) and /3|0 + ) — a|3~) (the higher energy barrier) with a > (3 > 0. 

Figure 10(c) corresponds to the case where the excitation energy is set equal to 
the barrier curvature, Hf2, which is 4.25 MeV in the present calculations. In this 
case, the lower energy barrier has an appreciable weight although the weight factor 
for the higher energy barrier is not negligible. When the excitation energy is further 
increased, the weight for the lower energy barrier becomes close to unity as is shown 
in Fig. 10(d), and the fusion cross sections are approximately given by 

a ins (E) = a fus (E; V(r) + A (r)), (4-12) 

where V(r) + An(r) is the lowest eigen-barrier (see Eq. (4-2)). Therefore, the main 
effect of the coupling to a state with a large excitation energy is to simply introduce 
an energy-independent shift of the potential, V(r) — > V(r) + Xo(r). This phenomenon 
is called the adiabatic potential renormalization. 96 ' 98 ) Typical examples in nuclear 
fusion include the couplings to the octupole vibration in 16 at 6.13 MeV 99/ and to 
giant resonances in general. 

In Refs. 97) and 98), it has been argued based on the path integral approach 
to multi-dimensional tunneling that the transition from the sudden tunneling to the 
adiabatic tunneling takes place at the excitation energy around the barrier curvature, 
hfl. That is, if the excitation energy is much larger than the barrier curvature, 
the channel coupling effect can be well expressed in terms of the adiabatic barrier 
renormalization. The numerical calculations shown in Fig. 10 are consistent with 
this criterion. 

§5. Fusion at deep subbarrier energies and dissipative tunneling 

Although the coudpled-channels approach has been successful for heavy-ion re- 
actions, many new challenges have been recognized in recent years. One of them is 
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Fig. 11. Fusion cross sections for the 64 Ni+ 64 Ni and 16 O+ 208 Pb systems as a function of the inci- 
dent energy. The experimental data are taken from Refs. 100) and 104). The dotted and the 
dashed lines are the results of potential model and the standard coupled-channels calculations, 
respectively. The solid lines denote the result when the deep subbarrier fusion hindrance is 
described in the adiabatic model. 108 ' 



the surface diffuseness anomaly discussed in Sec. 2.1. Another challenge, which may 
also be related to the surface diffuseness anomaly, 46 - 1 is an inhibition of fusion cross 
sections at deep subbarrier energies. This is a phenomenon found only recently, when 
fusion cross sections could become measured for several systems down to extremely 
low cross sections, up to the level of few nano barn (nb). 100 ^ 104 ^ These experimental 
data have shown that fusion cross sections systematically fall off much more steeply 
at deep-subbarrier energies with decreasing energy, compared to the expected energy 
dependence of cross sections around the Coulomb barrier. That is, the experimental 
fusion cross sections appear to be hindered at deep-subbarrier energies compared to 
the standard coupled-channels calculations that reproduce the experimental data at 
subbarrier energies, although the fusion cross sections are still enhanced with respect 
to a prediction of a single-channel potential model. 

Two different models have been proposed so far in order to account for the deep 
subbarrier fusion hindrance. In the first model, assuming the frozen densities in 
the overlapping region (i.e., the sudden approximation), Misicu and Esbensen have 
introduced a repulsive core to an internucleus potential, which is originated from 
the Pauli exclusion principle. 105 ^ See also Ref. 106) for a related publication. The 
resultant potential is much shallower than the standard potentials, and hinders the 
fusion probability for high partial waves. In the second model, on the other hand, 
Ichikawa, Hagino, and Iwamoto have proposed an adiabatic approach by assuming 
a neck formation between the colliding nuclei in the overlap region. 107 )> 108 ) i n this 
model, the reaction is assumed to take place slowly so that the density distribution 
has enough time to adjust to the optimized distribution. In this adiabatic model, 
the fusion hindrance originates from the tunneling of a thick one-body potential due 
to the neck formation. This model has achieved a comparable good reproduction of 
the experimental data to the sudden model, as is shown in Fig. 11. 
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The mechanism for the deep-subbarrier fusion hindrance has not yet been fully 
understood, as the two different models, in which the origin for the deep sub-barrier 
hindrance is considerably different from each other, account for the experimental data 
equally well. However, there is a certain thing which can be concluded by analyzing 
the threshold behavior in deep subbarrier fusion, 100 )' 101 )' 109 )~ m ) independent of the 
fusion models. no ) In Refs. 100), 101) and 109), the deep-subbarrier fusion hindrance 
has been analyzed using the astrophysical S factor. It has been claimed that deep- 
subbarrier fusion hindrance sets in at the energy at which the astrophysical S factor 
reaches its maximum. The authors of Refs. 100), 101) and 109) even parametrized 
the threshold energy as a function of charge and mass numbers of the projectile and 
the target nuclei. The relation between the threshold for deep-subbarrier hindrance 
and the maximum of the S factor is not clear physically, and thus it is not trivial 
how to justify the identification of the threshold energy with the maximum of the 
astrophysical S factor. Nevertheless, it has turned out that the threshold energy 
so obtained closely follows the values of phenomenological internucleus potentials at 
the touching configuration. no ) This strongly suggests that the dynamics which takes 
place after the colliding nuclei touch each other somehow makes the astrophysical 
S factor decrease as the incident energy decreases, leading to the fusion hindrance 
phenomenon. Notice that the fusion potential is almost the same between the sudden 
model and the adiabatic model before the touching (see Fig. 1 in Ref. 110)). 

One important aspect of fusion reactions at deep sub-barrier energies is that the 
inner turning point of the potential may be located far inside the touching point 
of the colliding nuclei (see Fig. 1). After the two nuclei touch each other, many 
non-collective excitations of the unified one-body system are activated. As is well 
known from the Caldeira-Leggett model, couplings to those excitations lead to en- 
ergy dissipation, which inhibits the tunneling probability. 10 ) The energy dissipation 
may occur also before the touching as a consequence of particle transfer processes to 
highly excited states in the target nucleus. 112 ) The phenomenon of deep subbarrier 
fusion hindrance may therefore be a realization of dissipative quantum tunneling, 
which has been extensively studied in many fields of physics and chemistry. A char- 
acteristic feature in nuclear fusion, which is absent or may not be important in 
dissipative tunneling in other fields, is that the couplings to (internal) environmen- 
tal degrees of freedom gradually set in. 121 ) That is, before the touching the fully 
quantum mechanical coupled-channels approach with couplings to a few collective 
states of separate nuclei is adequate, which however gradually loses its validity after 
the touching point due to the dissipative couplings. 104 ) This is the region in which 
the conventional coupled-channels approach does not treat explicitly by introducing 
an absorbing potential or by imposing the incoming wave boundary condition. Al- 
though it is highly important to construct a model for nuclear fusion by taking into 
account the dissipative couplings 108 )' 113 )' 114 ) in order to clarify the deep subbarrier 
fusion hindrance, it is still a challenging open problem. To this end, a transition from 
the excitations of two separate nuclei in the entrance channel, which are included 
in the conventional coupled-channels calculations, to molecular excitations (i.e., the 
excitations of the combined mono-nuclear system) has to be described in a consistent 
and smooth manner. 108 )' 115 )~ 117 ) The development of quantum mechanical version 
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of phenomenological classical models for deep inelastic collisions (DIC), such as the 
wall and window formulas for nuclear friction, 118 ^ 120 ) will also be important in this 
respect. 

§6. Application of barrier distribution method to surface physics 

The barrier distribution method discussed in Sec. 4 is applicable not only to 
heavy-ion subbarrier fusion reactions but also to any multi-channel tunneling prob- 
lem. In general, the barrier distribution is defined as the first derivative of penetra- 
bility with respect to energy, dP/dE (see Eq. (4-7)). 

As an application of the barrier distribution method developed in nuclear physics 
to other fields, let us discuss a dissociative adsorption process of diatomic molecules 
on a metal surface. When molecular beams are injected on a certain metal, such as 
Cu and Pd, diatomic molecules are broken up in the vicinity of metal surface to two 
atoms due to the molecule-metal interactions before they stick to the metal. This 
process is referred to as dissociative adsorption, and has been extensively studied in 
surface science together with the inverse process, that is, associative desorption. 122 ) 
The adsorption process takes place by quantum tunneling at low incident energies, 
as there is a potential barrier between the two phases of the molecules, i.e., the 
molecular phase and the breakup phase with two separate atoms. 122 )' 123 ) The vi- 
brational and rotational excitations of diatomic molecules play an important role 
in dissociative adsorption, 124 ^ 126 ) as in heavy-ion subbarrier fusion reactions. The 
coupled-channels method has been utilized to discuss the effects of the internal ex- 
citations of molecules on dissociative adsorption. 127 ^ 134 ) 

In this section, we discuss only the simplest case, that is, the effect of the rota- 
tional excitation on dissociative adsorption, while the vibrational degrees of freedom 
is assumed to be frozen in the ground state. In contrast to heavy-ion fusion reactions, 
the initial rotational state in the problem of dissociative adsorption is not necessarily 
the ground state. The initial rotational state of diatomic molecules can be in fact 
selected in molecular beams, and the experimental data of Michelsen et a/. 125 )' 126 ) 
have indicated that the adsorption probability of D2 molecules on Cu surface shows 
a nonmonotonic behavior as a function of the initial rotational state. That is, at a 
given incident energy, starting from the initial rotational state Lj = 0, the adsorption 
probability first decreases for = 5 and then increases for = 10 and Lj = 14 (see 
Fig. 9 in Ref. 126)). 

In order to explain this behavior, Difio, Kasai, and Okiji have considered a simple 
Hamiltonian for H2 and D2 molecules given 

by 130),132) 

where s is the one-dimensional reaction path in the two-dimensional potential energy 
surface spanned by the molecule-surface distance, Z, and the interatomic distance, r. 
The reaction takes place from s = —00, that corresponds to the approaching phase 
of molecules, to s = +00, where the incident molecule has broken up to two atoms. 
The sticking probability to the metal surface is identified as the penetrability of the 
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Fig. 12. The potential energy for the dissociative adsorption process of H2 molecule on metal surface 
given by Eq. (6-4). The parameters are £ a =0.536 eV, Vi = 1.0 eV, a=1.5 A~\ fi = 0.25, 
ro=0.739 A, and / = 0.14. The left panel shows the potential for L = as a function of the 
reaction path coordinate s for 6 = (the dashed line) and 9 = 7r/2 (the solid line), where 
6 is the molecular orientation angle (6=0 and 6 — tv/2 correspond to the configurations with 
perpendicular and parallel to the metal surface, respectively) and L is the associated angular 
momentum operator. The spherical part of the potential, Vo(s), is also shown by the dotted 
line. The right panel shows the sum of spherical part of the potential, Vq(s), and the rotational 
energy, H TOt (s) = L(L + l)ft 2 /2/(s), for three different values of L. 



potential barrier, V. M in Eq. (6-1) is the mass for the translational motion of the 
diatomic molecule given by M = 2m, where m is the mass of the atom (i.e., m = mu 
for H2 molecule and m = mo for D2 molecule). 6 is the molecular orientation angle, 
where 9=0 corresponds to the configuration of the molecule perpendicular to the 
surface while 9 = tt/2 to the configuration parallel to the surface. L is the associated 
angular momentum operator. J(s) is the momentum inertia for the rotational motion 
given by 

I(s) = f xr 2 (l + fe as ), (6-2) 

where [i = m/2 and / is a parameter characterizing the s dependence of the inter- 
atomic distance r, ro being the interatomic distance for an isolate molecule. The 
same parameter a as in Eq. (6-2) appears also in the potential energy, V(s, 9), which 
is parametrised as 

V(s,9) = ^ — -(l-/3cosfl 2 ) + yicos 2 6<- i(l + tanh(as)), (6-3) 

cosh (as) 2 

= V (s) + V 2 (s)Y 20 (9), (6-4) 



with 



V o(s) = 2 ° - ( 1 - £ ) + -£(1 + tanh(as)), (6-5) 



cosh 2 (as) V 3 J 6 

E a 2/3 [4tt 1 /4tt 1 
cosh 2 (as) ' T V T + 3 V T 



y 2(s) = ~ — 2 ° ~ ■ -r\l — + o \l —^i 1 + tanh(as)). (6-6) 



The coupled-channels equations for the Hamiltonian (6-1) can be derived in 
the same manner as in Sec. 3. For scattering with the initial rotational angular 
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momentum of molecules of Lj and its z-component Mj, we expand the total wave 
function as _ 

L 

Notice that the Hamiltonian (6-1) conserve the value of Mj, as the coupling potential 
is proportional to l2o(#)- The coupled-channels equations then read, 



ft2 * ■ £ < £+1 >*. +w .)- B 



2Mc2s 2 27( S ; 
where the matrix element ( Ylm, | Y20 \ Yy Mi ) is given by 
(Y LMl \Y 2 o\Y LIMi ) = (-) MI /JV(2L + 1)(2L' + 1) ( 



fe(5)+y2(s)^(yLM i |r 20 |^M i )fe i ( s ), 

(6-8) 



L 2 L' \ f L 2 L' 
) \ -Mi Mi 

(6 ' 9) 

Noticing that I(s) — > /itq for s — > — 00 and J(s) — > for s — > 00, these coupled- 
channels equations are solved by imposing the boundary conditions of 



<PllM = e*»8 L , Li - ^R LLi e~^ s (s -00), 



(6-10) 



^T LLi e lks (*->oo), (6-11) 

where k L = y / 2M(E - e L )/h 2 with e L = L(L + l)h 2 /2fir% and k = ^2MEjW . The 
adsorption probability for a given value of Lj and Mj is then obtained as 

PL iMi = Y,\ T ^f- ( 642 ) 

L 

By making average over all possible Mj, the total adsorption probability for Lj is 
given by 

^ = 2l^lE^" (^3) 

Let us now solve the coupled-channels equations for H2 molecules. The results 
are qualitatively the same also for D2 molecules. Following Ref. 130), we take 
£ a =0.536 eV, Vi=1.0 eV, a=1.5 A" 1 , /3 = 0.25, and r =0.739 A. For the factor / 
in Eq. (6-2), we take / = 0.14. 135 ) The potential with these parameters are shown 
in Fig. 12. The left panel shows the potential energy V(s, 9) given by Eq. (6-4) for 
two different values of 6. For comparison, the figure also shows the spherical part 
of the potential, Vq(s). One can see that the barrier is lower for the configuration 
parallel to the metal surface (that is, 9 = tt/2) as compared to the configuration 
perpendicular to the surface, 9 = 0. The right panel, on the other hand, shows the 
sum of the spherical part of the potential, Vq(s), and the rotational energy, H rot (s) = 
L(L + l)h 2 /2I(s), for three different values of L. Because of the s dependence of the 
rotational moment of inertia, I(s), the barrier height for the molecules incident from 
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Fig. 13. The results ol the coupled-channels calculation for the dissociative adsorption process of 
H2 molecules. The left panel shows the adsorption probability, P, while the right panels shows 
the barrier distribution defined as dP/dE for several values of the initial angular momenta Li 
for the rotational state of the molecule as a function of the initial kinetic energy -Ekm- 




Fig. 14. The results of the single-channel calculation, obtained by turning off all the coupling 
matrix elements in the coupled-channels equations except for the diagonal component. The left 
and the right panels show the adsorption probability, P, and the barrier distribution, dP/dE, 
respectively. The dotted lines denote the results when the initial rotational state is at Li = 0. 
The thin solid, the dot-dashed, and the thick solid lines are the results of (Li, M%) = (2, 0), (2, 1), 
and (2,2), respectively. The dashed lines with the solid circles show the results for Li — 2 
obtained by averaging all the Mi components. For the barrier distributions shown in the right 
panel, the weight factors, 1/5 (for M; = 0) and 2/5 (for = 1 and 2), are multiplied. 



s = —00, that is, the difference between the energy at s = and that at s = —00 
decreases as a function of L. 

The results of the coupled-channels calculations are shown in Fig. 13 for several 
values of initial rotational state, Li, in which the adsorption probability is plotted as 
a function of the incident kinetic energy of the molecule, defined as E = E^\ a +Li(Li+ 
l)h 2 /2rl. As has been noted in Refs. 130) and 132), these calculations well account 
for the non-monotonic behavior of the adsorption probability as a function of Lj. The 
right panel shows the corresponding barrier distribution, dP/dE, obtained with the 
point difference formula with the energy step of 0.03 eV. One can clearly see different 
structures for each Li. For Li = 0, the barrier distribution has three prominent peaks. 
These peaks are smeared for Lj = 4, and at the same time, the center of mass of 
the distribution is shifted towards high energy, leading to the decrease of adsorption 
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probability. This is due to the fact that the result for Lj = 4 is actually given by the 
average over contributions from 9 different Mj values. In order to demonstrate this 
effect, Fig. 14 shows the results of the single-channel calculations for Lj = 2 with 
three different values of Mj and its average. For comparison, the figure also shows the 
single-channel calculation for Lj = (the dotted line). We define the single-channel 
calculation as the one which neglects all the coupling terms in the coupled-channels 
equations (6-8) except for the diagonal term, L = L' . Because of the properties of 
the spherical harmonics, the diagonal term of the coupling potential is attractive for 
Mi = 2, while it is repulsive for Mj = and 1 (see Eq. (6-9)). The single peak 
in the barrier distribution for Li = is then distributed to three peaks in the case 
of Li = 4, shifting the center of mass of the distribution slightly towards higher 
energy (notice that — Mi gives the same contribution as Mj). With the off-diagonal 
components of the coupling potential, the distribution will be further smeared, as in 
the distribution for Lj = 4 shown in Fig. 13. When the initial angular momentum 
is further increased, the barrier distribution starts moving towards lower energies, 
as seen in the figure for Lj = 10 and 14, which enhances the adsorption probability 
as its consequence. This is mainly due to the fact that the barrier is lowered for a 
large value of rotational state, Lj, as has been shown in Fig. 12. 

The barrier distribution representation of the tunneling probability provides 
a useful means to understand the underlying dynamics of dissociative adsorption 
process, as the shape of the distribution strongly reflects the molecular intrinsic 
motions. This would be even so, particularly when the rotational and the vibrational 
degrees are taken into account simultaneously. 131 )' 133 ** It will be an interesting future 
study to investigate how the barrier distribution behaves in the presence of the 
rotational excitation together with the vibrational excitation. 

§7. Summary and outlook 

Recent developments in experimental techniques have enabled high precision 
measurements of heavy-ion fusion cross sections. Such high precision experimen- 
tal data have elucidated the mechanism of subbarrier fusion reactions in terms of 
quantum tunneling of systems with many degrees of freedom. In particular, the ef- 
fects of the coupling of the relative motion between the target and projectile nuclei 
to their intrinsic excitations have been transparently clarified through the barrier 
distribution representation of fusion cross sections. 

The channel coupling effects can be taken into account most naturally with 
the coupled-channels method. When the excitation energy of an intrinsic motion 
to which the relative motion couples is zero, the concept of barrier distribution 
holds exactly. In this case, quantum tunneling takes place much faster than the 
intrinsic motion. The effects of the couplings can then be expressed in terms of the 
distribution of potential barriers, and the fusion cross sections are given as a weighted 
sum of fusion cross sections for the distributed barriers. The underlying structure 
of the barrier distribution can be most clearly investigated when the first derivative 
of barrier penetrability, dP/dE, is plotted as a function of energy. In heavy- ion 
fusion reactions, this quantity corresponds to the second derivative of Ea{ us , which 
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is referred to as fusion barrier distribution, Df us . The fusion barrier distribution has 
been extracted for many systems through the high precision experimental data of 
fusion cross sections, <7f us . 

Even when the excitation energy of the intrinsic motion is not zero, the concept 
of fusion barrier distribution can be approximately generalized, using the eigen- 
channel representation of nuclear S-matrix, defined as the eigen-states of S^S. We 
have demonstrated that the barrier distribution shows a transition from the sudden 
to the adiabatic tunneling limits in a natural way as the excitation energy increases, 
where the potential is simply renormalized in the latter limit without affecting the 
shape of barrier distribution (i.e., the adiabatic barrier renormalization) . 

The barrier distribution representation is applicable also to other multi-channel 
quantum tunneling problems. A good example is the dissociative adsorption phe- 
nomenon in surface science. The rotational and vibrational excitations of diatomic 
molecules play an important role in the adsorption process. These effects can be de- 
scribed by the coupled-channels approach, and the barrier distribution can be defined 
as in heavy-ion subbarrier fusion reactions. The results of coupled-channels calcu- 
lations have indicated that the barrier distribution representation provides a useful 
means in clarifying the underlying mechanism in the dynamics of surface interaction 
of molecules. 

Although our understanding of subbarrier fusion reactions has considerably in- 
creased in the past decades, there are still many open problems in heavy-ion subbar- 
rier fusion reactions. For example, it has not yet been understood completely how 
the hindrance of fusion cross sections with respect to the standard coupled-channels 
calculation takes place at deep subbarrier energies. A promising mechanism of the 
hindrance is that many non-collective channels are activated after the target and 
the projectile nuclei overlap with each other, and the relative energy is irreversibly 
dissipated to the intrinsic motions. This would occur only at deep subbarrier ener- 
gies, in which the inner turning point of the potential barrier is located inside the 
touching radius of the two nuclei. This phenomenon may thus be a good example of 
dissipative quantum tunneling, which has been extensively discussed in many fields 
of physics and chemistry. A unique feature in nuclear physics is that the dissipative 
nature of the couplings gradually sets in, in a sense that the coupling is reversible 
before the touching and it gradually reveals the irreversible character as the overlap 
of the colliding nuclei increases. In order to gain a deep insight of this problem, 
it might be helpful to revisit heavy-ion deep inelastic collisions (DIC) from a more 
quantum mechanical point of view. This is important also in connection with the 
synthesis of superheavy elements by heavy ion collisions with large mass numbers, 
for which fusion cross is strongly hindered at energies near the bare Coulomb barrier. 

Other important issues not covered in this paper include fusion of halo nuclei 
and the role of multi-nucleon transfer. For the former, there have been many debates 
concerning how the breakup process affects subbarrier fusion. 136 H 42 ) However the 
interplay between fusion and breakup involve many complex processes 24 ) and the 
role of breakup in fusion has not yet been understood completely. Moreover, particle 
transfer processes also affect both fusion and breakup in a non-trivial way, as has 
been found recently in 6 ' 7 Li + 208 Pb reactions 143 ) (see Ref. 144) for a review on sub- 
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barrier fusion of weakly bound stable nuclei, 6 ' Li and 9 Be). A theoretical calculation 
has to take into account the fusion, transfer, and breakup processes simultaneously 
in a consistent manner. It remains a challenging problem to carry out such calcu- 
lations, although the time-dependent wave packet approach 141 ) has been performed 
with a limited partition for the transfer channels. From experimental side, fusion 
cross sections for many neutron-rich nuclei do not appear to show any particular en- 
hancement or hindrance, 145 )~ 148 ) but recent experimental data for 12 > 13 > 14 > 15 C-|- 232 Th 
reactions have shown that the fusion cross sections are enhanced for the 15 C projec- 
tile as compared to those with the other C isotopes. 149 ) Again, several sort of transfer 
channels would have to be considered to understand the differences in the behavior 
of fusion cross sections. 143 )' 150 )~ 152 ) In particular, the multi-nucleon transfer process 
may play an important role in fusion of neutron-rich nuclei. Although there have 
been a few attempts to treat the multi-neutron transfer process in subbarrier fusion 
reactions, 153 )~ 157 ) it has still been a challenging problem to include in a full quantum 
mechanical manner the multi-nucleon transfer processes consistently with inelastic 
channels by taking into account also the final Q-value distribution of transfer. 

A much more challenging problem is to describe heavy-ion fusion reactions, thus 
many-particle tunneling, 158 ) from fully quantum many-body perspectives, starting 
from nucleon degrees of freedom. The time-dependent Hartree-Fock (TDHF) the- 
ory has been widely employed to microscopically describe nuclear dynamics. 159 )' 16 °) 
It has been well known, however, that the TDHF method has a serious drawback 
that it cannot describe a many-particle tunneling phenomenon. In order to cure this 
problem, Bonasera and Kondratyev have introduced the imaginary time propaga- 
tion. 161 )' 162 ) In this connection, we wish to mention that an alternative imaginary 
time approach, called the mean field tunneling theory, for quantum tunneling of 
systems with many degrees of freedom has been developed in Ref. 163). The mean 
field tunneling theory is a reformulation of the dynamical norm method for quan- 
tum tunneling, 86 )' 164 ^ which evaluates the non-adiabatic effect on the tunneling rate 
through the change of the norm of the wave function for the intrinsic space dur- 
ing the evolution along the imaginary time axis. The mean field tunneling theory 
has been applied to quantum mechanically discuss electron screening effects in low 
energy nuclear reactions, 163 ) while the dynamical norm method has been used to 
discuss the effects of nuclear oscillation on fission. 164 ) It would be an interesting 
challenge to develop a fully microscopic version of these methods and apply them to 
heavy- ion fusion reactions. More recently, Umar et al. have used the density con- 
strained TDHF (DC-TDHF) method to analyze heavy-ion fusion reactions. 165 )~ 167 ) 
Even though these microscopic approaches seem promising, these are based on cer- 
tain assumptions, such as a local collective potential with single-channel. It is thus 
not yet clear whether they are applicable to many-particle tunneling problems in 
general, such as two-proton radioactivity 168 )~ 173 ) an d alpha decays. 174 H 78 ) It would 
be an ultimate goal to develop a general microscopic theory which can describe sev- 
eral tunneling phenomena simultaneously, not only in nuclear physics but also in 
other fields of physics and chemistry. Such theory would naturally provide a way 
to describe the role of irreversibility (that is, the energy and angular momentum 
dissipations) as well as the density evolution after the touching in subbarrier fusion 
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reactions without any assumption for the adiabaticity of the fusion process. 
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Appendix A 

Relation between surface diffuseness and barrier parameters 

In this Appendix, we discuss the relation between the surface diffuseness pa- 
rameter a in a nuclear potential and the parameters which characterise the Coulomb 
barrer, that is, the curvature, the barrier height, and the barrier position. With such 
relation, one can estimate the value of a from empirical barrier parameters. 

For a given nuclear potential Viv(r), the barrier position R b is obtained from the 
condition that the first derivative of the total potential is zero at r = Rb, 

~dV N {r) Z P Z T e 2 ^ 



dr r 



2 



= 0. (A-l) 

lr=R b 



r=R b 

The barrier height Vb and the curvature Q are then evaluated as 

V b = V N (R b ) + ^f^, (A-2) 
Kb 

a = ^p2M±|^fM, (A . 3) 

where V^(r) is the second derivative of the nuclear potential with respect to r. 

A.l. Exponential potential 

We first consider an exponential potential given by 

V N {r) = V e~ r / a . (A-4) 

From Eq. (A-l), the depth of the nuclear potential, Vq, is related to the charge 
product ZpZx as 

^ = 0. (A.5) 
a R b 

From this equation, the barrier height and the curvature read 
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respectively. 

A. 2. Woods-Saxon potential 

We next consider a Woods-Saxon potential given by 

V N {r) = V \ . . . (A-8) 

Combining Eqs. (A-l),(A-2), and (A-3), one finds that the surface difiuseness pa- 
rameter a is expressed in terms of Rb, Vb and Q as 

a ~ ^Q 2 Hl 2Z P Z T e 2 9 ) 

Z P Z T e 2 Z ZpZTe 2 -^,,^ 

Once the surface diffuseness parameter is thus evaluated, the other two parameters 
in the nuclear potential can be obtained as 



x ^ 2 Z P Zre 2 
Rl 



V = ae-%l + e*) 2 ^f-, (A-ll) 



where x is defined as (Rb — Ro) /a. 

Appendix B 

Parabolic approximation and the Wong formula 

If the Coulomb barrier is approximated by a parabola, 

V(r) ~ V b - ^Q 2 (r - R b f, (B-l) 
the corresponding penetrability can be evaluated analytically as 

P( E ) = ft-, \T" ( B ' 2 ) 

l + exp[fa{V b -E)] 

Using the parabolic approximation, Wong has derived an analytic expression for 
fusion cross sections. 49 ) He assumed that (i) the curvature of the Coulomb barrier, 
hfi, is independent of the angular momentum /, and (ii) the position of the Coulomb 
barrier, Rb, is also independent of I, and the dependence of the penetrability on the 
angular momentum can be well approximated by shifting the incident energy as 
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Fig. 15. (The left panel) The Coulomb barrier for the 16 0+ 144 Sm system shown in Fig. 1 (the solid 
line) and its parabolic approximation (the dashed line). (The right panel) Comparison of the 
corresponding fusion cross sections obtained by numerically solving the Schrodinger equation 
without resorting to the parabolic approximation (the solid line) and those obtained with the 
Wong formula, Eq. (B-5). 



If many partial waves contribute to fusion cross section, the sum in Eq. (2-17) may 
be replaced by an integral, 

7T r°° 

m m (E) = p dl (21 + l)Pi(E). (B-4) 

Changing the variable from I to 1(1 + 1), the integral can be explicitly carried out, 
leading to the Wong formula 49 ) 



a fus (£) = — Rlln 



(B-5) 



At energies well above the Coulomb barrier, this formula reduces to the classical 
expression of fusion cross section given by Eq. (4-3). 

The left panel of Fig. 15 shows the parabolic approximation to the Coulomb 
barrier for the 16 O + 144 Sm system shown in Fig. 1. Because of the long ranged 
Coulomb interaction, the Coulomb barrier is asymmetric and the parabolic potential 
has a smaller width compared with the realistic potential. Nevertheless, the Wong 
formula for fusion cross sections, Eq. (B-5), works well except at energies well below 
the barrier, where the parabolic approximation breaks down (see the right panel of 
Fig. 15). 

Even though the Wong formula appears to work well for the single-channel 
potential model, one can still discuss the corrections to it. The first correction 
is with respect to the integral in Eq. (B-4). To discuss the correction, we first notice 
that replacing the sum in Eq. (2-17) with the integral in Eq. (B-4) is equivalent to 
taking only the leading term (m = 0) of the exact Poisson sum formula, 

n 00 roc 

aUE) = y I Y,( 2l + 1 ) P l( E )=tf E / AP(£;A)e 2 ™ iA dA, (B-6) 

l m=-oo ^° 

where P(E\ A) is any smooth function of A satisfying P(E, 1 + 1/2) = Pi(E). 81 ^ Poffe, 
Rowley, and Lindsay have evaluated the contribution of the next most important 
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terms, m = ±1. 179 ) These terms lead to an oscillatory contribution to the fusion 
cross sections, 

o-fu S (E) = a w {E) + a osc (E), (B-7) 
where aw(E) is given by Eq. (B-5), while the oscillatory part a osc (E) is given by 

a osc (E) = ^fiRl^ exp I f^T" I sin(27rZ 9 ). (B-8) 

Here, l g is the grazing angular momentum satisfying 

£=nr)+ «k±i)i (B . 9) 

For heavy systems, the oscillatory part of fusion cross sections, cr osc , is usually much 
smaller than the leading term, aw- However, for light symmetric systems such as 
12 C+ 12 C, the oscillatory part becomes significant. 179 )~ 183 ) For a system of identical 
spin-zero bosons, the factor (1 + (—1)') has to be multiplied in the sum in Eq. (2-17) 
due to the symmetrization effect, making the contributions of all the odd partial 
waves vanish. In this case, the leading term of the fusion cross section is still given 
by the Wong formula, Eq. (B-5), while the oscillatory part becomes 179 -* 

<r osc (E) = AkhRI™ exp (~^0^-) sin(7r/ fl ) . (B-10) 

Fig. 16 shows the fusion cross sections for the 12 C+ 12 C reaction obtained with a 
parabolic potential with Vb = 5.6 MeV, i4=6.3 fm, and hfi = 3 MeV. The solid 
line shows the result of the exact summation of partial wave contributions with Eq. 
(B-3), while the dashed line shows a sum of Eqs. (B-5) and (B-10). A separate 
contribution from the Wong formula, Eq. (B-5), is also shown by the dotted line. It 
is seen that the oscillation of fusion cross sections can be well reproduced with the 
formula given by Eq. (B-10). 

The second correction to the Wong formula is the angular momentum depen- 
dence of the barrier radius. 184 ) Up to the first order of h 2 //j, 2 f2 2 R'l ) , Balantekin et 
al. have shown that the barrier radius for the Z-th partial wave Ru is given by 

This equation indicates that the barrier position decreases as the angular momentum 
I increases. At energies well above the barrier, the classical fusion cross sections are 
then modified to 184 ) 

aUE) = nRl (l " I) " ^(E - V b ) 2 (E » V b ). (B-12) 

Comparison between Eqs. (4-3) and (B-12) shows that the Wong formula slightly 
overestimates fusion cross sections at energies well above the Coulomb barrier. 
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Fig. 16. The fusion excitation function for the C+ C system. The solid line is obtained by 
carrying out exactly the angular momentum sum (with the symmetrization factor) in Eq. (2-17) 
with a parabolic potential with Vb = 5.6 MeV, Rb~6-3 fm, and hfi = 3 MeV. The barrier position 
and the curvature are assumed to be independent of the angular momentum I. The dotted line 
is obtained with the Wong formula, Eq. (B-5), while the dashed line is obtained as a sum of the 
Wong formula and the oscillatory cross sections given by Eq. (BTO). The experimental data 
are taken from Ref. 182). 



Appendix C 

Multiphonon coupling 

In this Appendix, we show that the dimension of coupled-channels equations can 
be reduced for vibrational couplings by introducing effective multi-phonon channels. 
Suppose that we have two modes of vibrational excitations (e.g., a quadrupole and 
an octupole modes), and consider the excitation operator 

6 = p 1 {a\ + a 1 ) + fo{a\ + a 2 ), (C-l) 
and the phonon Hamiltonian 

Hq = h/jj\a\a\ + hjjj 2 a\a 2 , (C-2) 

where a\ and a} 2 are the phonon creation operator for the first and the second modes, 
respectively. ft (i = 1,2) are the coupling constants, while hwi (i = 1,2) are the 
phonon excitation energy for each mode. We have shifted the phonon energies so 
that the ground state is at zero energy. 

If we truncate the phonon space up to the one phonon states, we have three 
basis states, 1 00) , |10), and 1 01) , where the state \n1n2) corresponds to the product 
state of n\ phonon state for the first mode and n 2 phonon state for the second mode. 
Here we have included the states with n\ + n 2 < 1. The matrix elements of the 
operator Hq + O with these basis states read, 






Pi 


ih 


Pi 







(h 





fVjJ 2 



H + O = ft foji , (C-3) 
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It is easy to see that the ground state 1 00) couples only to a particular combination 
of 1 10) and 1 01), 69 ) 

|1)= J_ (A|10)+ft|01)), (C-4) 
V Pi + "2 

with 



O|00) = ^/3 1 2 + /3 2|i). (C-5) 

The other combination of 1 10) and |01), ^2! 10) — /?i 1 01) , couples neither to 1 00) nor 
|1), and this can be removed from the coupled-channels calculation if the excitation 
energies of the two modes are the same, hooi = fiW2 = In this case, the dimension 
of the coupled-channels equations can be reduced to two with a modified strength 

H « + 6 ={1 L)- < C ' 6 » 



where j3 is defined by /3 = \J 0\ + One can easily generalize this scheme to higher 
members of phonon states. The resultant matrix is equivalent to that for a single 
phonon mode with the effective strength f3. For instance, when the phonon space is 
truncated at the two-phonon states, the coupling matrix reads 




where the effective two phonon state is defined as 

|2> = ^2^2 (^l 20 ) + >/2ft&|ll> + /3 2 2 |02)) . (C-8) 
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